Finite-size scaling behavior in trapped systems 
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Numerical transfer-matrix methods are applied to two-dimensional Ising spin systems, in presence 
of a confining magnetic field which varies with distance \x\ to a "trap center", proportionally to 
(|x|/^) p , p > 0. On a strip geometry, the competition between the "trap size" I and the strip width, 
L, is analysed in the context of a generalized finite-size scaling ansatz. In the low-field regime I 3> L, 
we use conformal-invariance concepts in conjunction with a linear-response approach to derive the 
appropriate (p-dependent) limit of the theory, which agrees very well with numerical results for 
magnetization profiles. For high fields I < L, correlation-length scaling data broadly confirms an 
existing picture of p-dependent characteristic exponents. Standard spin-1/2 and spin-1 Ising systems 
are considered, as well as the Blume-Capel model. 
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I. INTRODUCTION 

The continuingprogress in the field of trapped atomic 
quantum gases [J Q has led to the development of a 
variety of experimental techniques (e.g., in situ imag- 
ing [3j) to investigate Mott insulating phases in trapped 
bosonic or fermionic atoms Q ■ Inspired by the physics of 
strongly-correlated systems, a further step in the agenda 
will be the investigation of spin-ordered phases, given the 
recent improvements on cooling methods and ultra-low- 
temperature measurements. An important ingredient of 
the latter advance is the nature and shape of the atomic 
density profile within the trap, leading to a free gas at the 
wings of the distribution 0, Q . As the accuracy of these 
measurements improves, a detailed description of the 
possible crossovers involved will be needed in order to an- 
alyze the experimental data. Indeed, for an atomic sam- 
ple of finite size, L, the presence of a trap brings about 
another length scale, I (to be defined more precisely be- 
low), characterizing the trap shape and size. Therefore, 
the competition between these (geometric) lengths and 
the one describing the decay of correlations, £, which is 
assumed to diverge at the ordering phase transition in 
the thermodynamic limit, must be incorporated in a new 
scaling theory, which generalizes the standard finite-size 
scaling (FSS) 0. We should mention that a trap-size 
scaling regime, in which both £ and I are large, has been 
considered in Ref. with the relevant parameters ad- 
justed in such a way that £ < L, thus eliminating most 
finite sample-size effects. 

Our purpose here is to develop such generalized scaling 
theory, and we also adopt the strategy of leaving aside all 
quantum effects Q and concentrating instead on classical 
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spins localized on sites of a lattice, subject to a trapping 
field. More specifically, we consider Ising-like spins un- 
der a spatially varying magnetic field, whose intensity in- 
creases with distance |a?| to the origin (trap center) [cllicj. 
This way, for suitably large \x\ the spins are pinned par- 
allel to the local field. With the well-known mapping of 
the (pseudo)-spin variable onto a lattice-gas picture, one 
thus gets at least the main qualitative features of trap- 
ping phenomena. 

We investigate the scaling behavior of two-dimensional 
Ising spins close to the critical temperature T c at which, 
in zero field and in the thermodynamic limit, the cor- 
relation length £ becomes much larger than the lattice 
parameter, a (to be taken as unity in what follows). Ini- 
tially, only spin— 1/2 systems are considered. Upon addi- 
tion of a trapping magnetic field, as outlined above, £ is 
prevented from diverging, thus smearing out the second- 
order phase transition. Nonetheless, we can still look for 
signatures of the smeared phase transition in both the 
trap scaling regime t < L, considered in Ref. as well 
as in a "shallow trap" regime, in which I ^> L. 

In both regimes we use numerical transfer matrix (TM) 
methods and FSS ideas to disentangle the respective 
trapping and finite-size effects in Ising-like systems which 
are necessarily finite in extent. The quantities of in- 
terest, such as free energies, equilibrium site magneti- 
zations, bond energies, correlation functions and corre- 
lation lengths, can be extracted from the leading eigen- 
values of the TM, and from their respective eigenvec- 
tors [ll|. As we will see, in the shallow trap regime, 
conformal invariance allows us to develop a perturbative 
theory, whose predictions are in excellent agreement with 
the numerical data. We investigate the quantitative de- 
pendence of scaling behavior on trap shape, and test for 
universality against spin magnitude, the latter by consid- 
ering both a standard spin-1 Ising system and the Blume- 
Capel model. 

The layout of the paper is as follows. In Sec. [TT| we 
parametrize the trapping field in more detail, and present 
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the scaling ansatz for the free energy, making contact 
with previous work Q, when applicable; highlights on 
some technical aspects of the methods used in this pa- 
per are also given. Section IIHI deals with the low-field 
(or shallow trap) regime, in which scaling predictions 
are drawn perturbatively from conformal invariance, and 
checked against numerical TM data. TM methods are 
used to probe the high-field (or steep trap) regime in 
Sees. IIVI and [V] the latter concerning the behavior of 
S = 1 spins. Finally, in Sec. I VII concluding remarks are 
made. 



II. BASIC ASPECTS 

We consider a trapping field with a single power-law 
dependence on distance to the trap center: 



h(x) = h \x\P = (\x\/£f 



(1) 



Thus £ = h Q is an effective "trap size". In cold-atom 
experiments, harmonic potentials {p = 2) are usual. 

We use TM methods on long strips of width L sites 
of a square lattice, with periodic boundary conditions 
(PBC) across. So £ 3> L corresponds, for all x inside 
the strip, to a low-field regime where at T near T c and 
large L the system is still in the scaling regime (Sec- 
tion [TITl); an d the same is true at |x| near zero, even for 
£ = O(L) ("high fields", Section HV). In the strip geom- 
etry, the trapping field is translationally invariant along 
the "infinite" direction; given Eq. (fTJ), consistency with 
PBC demands that the field be symmetric relative to 
the line halfway along the strip width, thus in this case 
the distance |x| of Eq. (fTJ) stands for position relative to 
such an axis. It is expected [1(1 [HI, and has been ver- 
ified numerically @, that the scaling properties of trap- 
ping in this one-dimensional well will be the same as in 
a full two-dimensional one, for which the equipotentials 
are concentric circumferences. 

With t = (T — T c )/T c , and a uniform magnetic field H, 
a FSS expression for the singular part of the free energy 
can be adapted to the trapping context, as follows: 



F(t,H,i,L) = b~ d [ F(tb Vt ,Hb VH \t~ x b yi ,lrH Vh ) 



where b is an arbitrary rescaling parameter; y t = 
!Jh = (d + 2 — ri)/2 are the usual scaling exponents; y^ = 
1, and the "trapping" exponent ye is given, in this case 
where the trapping field couples to the magnetization, 
by 1: 



2p 



(3) 



— Q 

ye~ d + 2-r) + 2p' 

Making b = L, one gets: 

F(t, H, L,£) = L- d F(Lt",HL VH ,L( 

This prescribes the scaling behavior, in the trap, of all 
thermodynamic properties. Beyond scaling, their spe- 
cific forms are not provided by Eq. ((4]), and they will be 
modified by the trap. 
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III. LOW-FIELD REGIME 

The trapped system we consider is a spin- 1/2 Ising 
one, (spin variables a = ±1 on sites of a square lat- 
tice) with nearest-neighbor interactions of unit value, 
so the exact critical point is at T c = 2/ln(l + s/2), 
H = l~ l = L- 1 = 0. Setting t = H = 0, we first fo- 
cus on the low- field regime corresponding to I ^> L. In 
renormalization-group terminology, one has two relevant 
nonzero fields, namely L _1 and As the effects of 

the first of these are well-understood via conformal in- 
variance [llj], they can be readily incorporated into our 
description. In the regime under consideration, the re- 
sulting picture is a suitable starting point for a pertur- 
bative treatment of the second field. 

Within a linear-response context, and with aid of the 
fluctuation-dissipation theorem, one can write for the 
equilibrium magnetization at site i: 
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(<Ti(Tj) c hj 



(5) 



where (• ■ • ) denotes thermodynamic average, hj repre- 
sents the trapping field at site j, and the subscript c 
stands for connected spin-spin correlation functions. 

In keeping to the first-order perturbation picture 
adopted here, one substitutes the zero- field expressions, 
Cij = (<7i (jj)° c on the right-hand side of Eq. (El). These, 
in turn, are given via conformal invariance [13l | , on a strip 
of width L with PBC, by: 



C 



(27T/L)" 



[2cosh(27rx/L) - 2cos(2 7 r?//L)] ?,/2 ' 



V = 1/4 ■ 



(6) 

Here, x, y are the relative spin-spin coordinates, respec- 
tively along the strip, and across it. 

Though strictly speaking, Eq. ^ is an asymptotic 
form, discrepancies are already very small at short dis- 
tances, amounting to less than 0.5% for separations of 
the order of three lattice spacings flil . [l5| . 

Setting the origin at a point on the trapping axis, thus 
hj = h(x,y) = ho \y\ p , using Eq. © in Eq. (JSJ, and 



(2) transforming the latter into an integral, one gets: 



Hv = o)) 



'I 



2+p-r/ 



A( P ) =A 1 (p) + A 2 ( P ) 



A( P ) , (7) 



(8) 



where A\ (p) represents the long-distance contribution 
(x/L > 1, say) for which it is justified to ignore the y— de- 
pendence of Eq. ((6]), and A 2 (p) takes fully into account 
the short-range, angle-dependent, part of the integral. 
For x, y <C L in the latter, the denominator in Eq. ^ is 
proportional to R 2 = x 2 + y 2 , and one can use polar co- 
ordinates for the integration close to the origin. Further 
contributions of order 1/L to A{p) have been neglected. 
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One finds: 



=W2 



-P+l 



rjix 



1 faTr)^^ 2 

it p — 1] + 2 



(9) 
(10) 



where a is a constant of order unity (reflecting the choice 
of region where polar coordinates are used) , and C (p) = 

£ /2 (sm9)P d6 = [ (p - 1)!! / pU] a, where a = 1 for odd 
p, and 7r/2 for even p. 

The leading scaling behavior of the magnetization can 
be extracted from Eq. by first replacing I = h Q 1//p 
there, thus showing that ho and L always occur in the 
combination ho L p / e in the free energy. However, as can 
be seen from Eq. (fTJ), the field amplitude ho does not have 
the dimension of a magnetic field. Accordingly, given 
the non- uniform character of h(r), the appropriate field 
variable (with the correct dimension) to use when dif- 
ferentiating the free energy to obtain the magnetization 
is a spatially averaged, coarse grained, one; for the strip 
geometry considered here it reads, 



L 



L/2 



L/2 



h \y\Pdy~h L*> 



(11) 



Finally, within a linear response approach, a linear de- 
pendence of the magnetization on field corresponds to a 
quadratic dependence of the free energy. We then have 
F ~ L~ d {hoLv/ ) 2 , which yields 



WT^dh hoL 



(12) 



in agreement with Eq. (Q. Note further that the above 
derivation is valid also away from the central site (y = 0); 
see below. 

The corresponding predictions for (ero) have been 
checked against TM results. In Fig. Q] we show (ctq) 
against strip width L (10 < L < 24), for p = 1 and 
2, and ho = 10 -5 in Eq. (fT]), which amounts to having 
£ = 10 5 and 316.2 respectively for p = 1 and 2. Sin- 
gle power-law fits to TM data give p + 2 — r/ = 2.752(9) 
(p = 1), and 3.70(3) (p = 2), illustrating that the agree- 
ment with Eq. ([7]) indeed improves as the ratio L/l — > 0. 

Furthermore, one can test the amplitude ratio 
A(pi)/A(p 2 ), by plotting (<7 (pi, L))/(a (p 2 , L)) against 
L Pl ~ P2 . Our results for p\ = 1, p 2 = 2 are shown in 
Fig. O TM data are well-adjusted by a linear least- 
squares fit (y = ax + b) with b = 0.001(1), providing 
complementary evidence in favor of the single power-law 
scenario of Eq. (JT]). The adjusted slope is a = 3.01(2), 
to be compared to the predictions of Eqs. © and (fTU|). 
A(l)/A(2) = 2.98(1), where the uncertainty follows from 
allowing a of Eq. (fTQ| to vary between 1.0 and 1.5. 

We conclude that the linear response approach, in con- 
junction with conformal invariance concepts, is a suitable 
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Figure 1: (Color online) Double-logarithmic plot of magne- 
tizations at strip center for p = 1 and p = 2, against strip 
width L. Here, ho = 10~ 5 . The straight lines are guides to 
the eye, with respective slopes 2.75 (p = 1) and 3.75 (p — 2); 
see Eq. 0. 
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Figure 2: (Color online) Plot of the ratio of magnetizations 
on the trapping axis, {ao(pi, L)} / (ao(p2, L)} against L P1 ~ P2 
for 10 < L < 24, p! = 1, p 2 = 2, with ho = 10~ 5 . Points are 
TM results; the line is a linear least-squares fit of data (see 
text). 
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Figure 3: (Color online) Scaled magnetization TM data, 
Hq 1 L- {p+2 - v) (a yi ), for L = 24 (points), together with the 
result of linear response plus conformal-invariance approach 
(solid line). Both for p = 1. The horizontal axis is u — yi/L. 



description of the low-field magnetization properties of 
the trapping problem on a strip geometry. 

Our next step is to investigate the magnetization pro- 
files across the system. If the origin of coordinates in the 
integrals is now set at a distance y\ from the trapping 
axis, < y\ < L/2, and defining the reduced coordinate 
u = yi/L, one can see that the long-distance contribu- 
tion Ai(u,p) to A(u,p) [= A{u), for short] is still given 
by Eq. The short-distance part ^(u) can no longer 
be analytically expressed in a relatively simple form, be- 
cause the symmetries used in establishing Eq. (fTU)) no 
longer hold off-axis. However, one can still evaluate the 
corresponding integrals numerically, in which case the 
full form of Eq. ^ can be used over the whole domain 
of integration. We have found that, in general: 

(i) the order of magnitude of Ai (u) remains similar to 
that of Ai, as is the case on the strip axis [see Eqs. (0) 
and (HSJ)]; 

(ii) the fractional variation (A(l/2)-A(0))/A(0)) from 
axis to edge increases with p, from ss 2.5% for p = 1 to 
f« 6% for p = 8; 

(iii) as suggested by Eq. ([7]), together with generic 
scaling ideas, one gets good scaling plots of TM results, 
l -(p+2- v ) a g a i ns t u . 

We first illustrate the suitability of the linear-response 
plus conformal-invariance approach, in the description of 
numerical TM data for magnetization profiles. Figure [3] 
shows scaled magnetization data for p = 1, L = 24, as 
well as A(u) calculated by numerical integration. Apart 
from an overall normalization factor for A, which affects 
only the vertical scale, there are no fitting parameters. 
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Figure 4: (Color online) Calculated profiles A(u) for p = 
1, 2, and 8. Points are scaled magnetization TM data, 
hg 1 L- (p+2 -^(a yi }, for p = 1 and 2, both for L = 24. The 
horizontal axis is u = yi/L. All curves have been normalized 
to unity at u — 0. 

One sees that the agreement is very good. 

As regards point (ii), Fig. [5] provides a visual com- 
parison of calculated profiles A(u) for p = 1, 2, and 8, 
together with TM results for p = 1 and 2 (the latter, 
both for L = 24). In order to exhibit fractional vari- 
ations more clearly, all curves have been normalized to 
unity at u = 0. Note that, for p = 2, the agreement be- 
tween TM data and the respective calculated profile is as 
good as that for p = 1 . TM calculations for the low-field 
regime L/^< 1 with p = 8 have proved rather unwieldy, 
as the corresponding fields on strip sites become numeri- 
cally extremely small, thus we do not present a compar- 
ison between the calculated profile and TM data for this 
case. In this connection, one must recall that the phys- 
ically more interesting trap shapes correspond to p = 1 
and 2, since p large approaches an infinite well. 

For the data collapse mentioned in point (iii) , we found 
that strong even-odd oscillations are present, thus in the 
following we restrict ourselves to even lattice widths. We 
also consider only p = 1, which gives the smoothest re- 
sults, as expected from the fact that this setup corre- 
sponds to a shallow trap shape. Finally, corrections to 
scaling are present which are small but noticeable for 
small and intermediate strip widths. These have been 
accounted for in the usual way [l6[, by assuming 

h -i L - {p+ 2- v ) Ki) = f p {y x /L)+L- u g 9 {y 1 /L) . (13) 

Here, to > is the exponent associated to the leading 
irrelevant operator. We found that, with u = 1, and 
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Figure 5: (Color online) p = 1: data collapse for scaled mag- 
netization TM results, including corrections to scaling, i.e., 
h^ 1 L- 2 - 75 {a Vl )-L-"gi(yi/L) [see Eq. JS)], for 10 < L < 24 
(even widths only). The horizontal axis is u = yi/L. 



gi(y±/L) = c (constant), the best data collapse is ob- 
tained for c = -7.0(2) x 1CT 5 . This is shown in Fig. [5] 



IV. HIGH-FIELD REGIME 

In order to make contact with previous work Q, we 
have set p = 2, and generated data for strips with varying 
L and £, initially keeping constant the ratio £/L = 1/ s/2. 

The average energy per site e s can be calculated, for 
a system with nearest-neighbor interactions on a square- 
lattice geometry, as g s — &x ~l~ where e x , 6y are the 
energies of the bonds along each coordinate axis. On a 
strip of width L sites with PBC across, with the TM ad- 
vancing along x, one finds for the Ising spin-1/2 model at 
T c and in the absence of (trapping or uniform) field, that 
e x, y = e^+a Xty L~ 2 + b x ^ y L~ 4 +0(L~ 6 ). Sequences of e x , 
e y generated using 4 < L < 16 give a v /a x = b y /b x = —1 
to four significant digits (a y ,b y > 0). Both extrapolate 
to the exact value e° = \/2/2 to four significant digits, 
while the average between the two extrapolates agrees 
with that to six digits. 

Upon inclusion of the trapping field as described above, 
its value at sites near the strip edge becomes rather large. 
The site energies now vary with position across the strip. 
Since signatures of a scaling regime should be more ap- 
parent where deviatons of the relevant fields from their 
critical values are small, at least in a local sense, we cal- 
culate energies for the sites at, or nearest to, the strip 



center. A similar reasoning appears to have been fol- 
lowed in the Monte Carlo simulations of Ref. H. 

The resulting sequences of e x and e y exhibit even-odd 
oscillations. However, these are smoothed out when e s 
is considered instead. Fitting data for 10 < L < 23 
to e s = y/2 + a£- e (1 + bT 2 ) gives 9 = 0.514(10), in 
very good agreement with the scaling prediction p], 6 = 
16/31 = 0.51613... . 

Next, we consider local magnetizations mo, again on 
sites at, or nearest to, the strip center, for the same rea- 
sons invoked above in regard to site energies. The magne- 
tization sequences show significant even-odd oscillations, 
thus we ran separate fits for data subsets correspond- 
ing to even and odd L. Using 10 < L < 23, assuming 
m = cl-^ (1 + d£- 2 ) we get = 0.0646(8) from an av- 
erage between the final estimates for even- and odd- se- 
quences. This is again in very good agreement with the 
respective scaling prediction [8|], <fi = 2/31 = 0.0645 .... 

In analogy with Eq. (j4]), and considering the uniform 
field H = from the start, the correlation length on a 
strip is given by: 

H L (tJ)=L.f(Lt 7 Lr 9 ) (14) 

At t = 0, and applying standard FSS ideas, one expects: 



L 



7T £((),*) 



{const. 
x 



X < 1 

x > 1 



(15) 



We generated correlation-length data for 4 < L < 24 
and assorted 6 < I < 20, thus allowing Ljl to vary 
over a relatively broad extent. Recalling from conformal- 
invariance [l3[ that L/n £l(Q,0) = r\ = 1/4, we plot- 
ted L/tt £l(0,£) against Li~ 9 for tentative values of 0. 
searching for the best data collapse. This was found along 
a somewhat broad range, 0.49 < 9 < 0.52, which includes 
both the trap-size scaling prediction and the mean-field 
value #mf = 1/2- Fig. O shows a typical result, using 
6 = 0.51. The trends predicted in Eq. (JTSj are all veri- 
fied. 

In order to narrow down our error estimates, we inves- 
tigated how well the large-x portion of the data falls on 
a straight line, as predicted in Eq. dT3J - We fitted por- 
tions of data contained within intervals x\ < x < X2 to 
a single power law with an adjustable exponent, e. For 
[£1,22] = [4,9], i.e. essentially all the large- a; region, we 
found e = 1 within 0.5% for 9 = 0.503(7), while if we 
restricted ourselves to [x 1,^2] = [4,6] the corresponding 
range of 9 turned out to be 9 = 0.514(10). 

Since the largest values of x correspond to large £ _1 , 
the respective correlation-length data belong to a region 
in parameter space away from the critical point, where 
the scaling behavior crosses over to a mean-field picture. 
Accordingly, when one includes such data in the analysis, 
the results for 9 are biased towards its mean-field value; 
on the other hand, their exclusion brings the estimate of 
9 back towards the scaling prediction. 

An analytic treatment of the correlation length can 
be given for the limit p — > 00, providing for this case 
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Figure 6: (Color online) Scaled correlation-length data for 
p = 2, with varying L/£ [see Eq. (|15|l ]. The horizontal axis is 
x = with # = 0.51. The straight line through the origin 

is a guide to the eye. 

the scaling function g(x) in Eq. (|15j) . and having simi- 
larities to the form shown for p = 2 in Figure [H] The 
reduction for p — > oo comes because the trapping well is 
there equivalent to confinement between —£ and +£, so 
strip geometry and well are together equivalent to a strip 
of width L' = min(L, 2£). For this equivalent system, 
Eq. © (with L replaced by L') can be used to obtain 
the effective correlation length, yielding for p — > oo, 

TT^O.O = 5 ' - 9( ' T) = 2 maX(2 ' ^ (16) 

[using 6(p = oo) = 1, from Eq. ©]. 

V. UNIVERSALITY: ISING S = 1 AND 
BLUME-CAPEL MODEL 

In this section, we investigate the universality of trap- 
ping scaling against varying features of the spin Hamilto- 
nian. We shall keep p = 2 in what follows. For TM cal- 
culations, the field setup is the same as defined in Sec.HIl 
with the trapping axis halfway along the strip width. 

We start by applying the ideas discussed in the pre- 
vious Section to an Ising 5 = 1 system, on a square 
lattice with nearest-neighbor couplings, in units of which 
the zero-field critical temperature is T c = 1.69356(2) [l7| . 
With the trapping field coupled to the magnetization, we 
considered the high- field regime, making again £ = L/y/2 
as our starting point. 

At the zero- field critical point, by considering the bond 
energies e x and e y as in Sec. IIV1 we found from fits 



of data for 4 < L < 13 that e° = 1.16100(2), to be 
compared with e° = 1.16094(5) [lq . Upon introducing 
the trapping field, we restricted ourselves to the sites 
closest to the strip center, for the reasons invoked in 
Sec. [TV] Fitting data for 4 < L < 15 (even- and odd-L 
sequences have to be considered separately, on account of 
the associated oscillations) to e s = + ai£~ (l + b\£~ 2 ) 
gives 9 = 0.55(2), where the error bar reflects the 
spread among fits of differing subsets of data: odd/even, 
and vertical/horizontal/ vertical-plus- horizontal bond en- 
ergies. 

For the local magnetization closest to the strip cen- 
ter, for 4 < L < 15 we again resorted to separate fits 
for even- and odd-L sequences. The relative shortness of 
such sequences, compared to the 5 — 1/2 case, makes 
extrapolation somewhat risky. However, by attempting 
power-law fits too = ci£ _< *(l + di£~ 2 ). we saw that, by 
successively excluding an increasing number of smaller 
strips widths, the sequences of estimates for <f> vary as: 
(f> = 0.088(4), 0.075(2), 0.070 (L even); and = 0.099(6), 
0.078(2), 0.071(2), 0.068 (L odd). Thus, it seems plausi- 
ble that consideration of larger L would produce a result 
compatible with the prediction = 2/31 Q. 

Next, by allowing the trapping field intensity (thus, 
the corresponding length £) to vary independently of L, 
we extracted an independent estimate of 9 via the scal- 
ing argument presented in Eqs. (fT4"f and (|T5|) . We used 
4 < L < 15, and 2 < £ < 9. For tentative values of 6 
close to 0.5, we found a scaling curve quantitatively very 
similar to that for S = 1/2, shown in Fig. [5] Indeed, the 
respective numerical values stay within at most 3 — 4% 
from those of the 5 = 1/2 curve. However, the quality 
of data collapse is markedly inferior; furthermore, it does 
not appear possible to have good scaling with a single 6 
for the whole curve. For = 0.47, points fall very well on 
a straight line at large x = L/£ e , while those at smaller 
x show a somewhat large degree of scatter. For 6 = 1/2 
the situation is reversed. 

Overall, the above results are in fair agreement with 
the hypothesis of universal behavior against spin value, 
though with some degree of scatter. Note that the esti- 
mates for 9 from internal-energy and correlation-length 
scaling are, respectively, above and below the expected 
value 9 = 16/31; thus, we have found no systematic trend 
away from universality. 

We now turn to discussing the application of a trapping 
field to a variant of the 5 = 1 Ising model, namely the 
Blume-Capel (BC) model [l9l - l23| . on a two-dimensional 
square lattice. The Hamiltonian is: 

U = -jY J SiS ] + (D + h \n\P)Y^S 2 l , (17) 

where Si = ±1,0 (spin-1 Ising variables). With J > 
(taken to be unity from now on) and D, fig > 0, the 
field term privileges vacancies, so there is an immediate 
mapping between (5 2 ) and the "particle" density in the 
lattice-gas picture. The BC model was studied recently 
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via Monte Carlo simulations,_in a trapping field with p = 
1 , and assorted values of D . 

Considering for now ho = 0, it is known that for D > 
the competition between field and spin-spin coupling 
leads to a tricritical point at T t — 0.6085776(1), D t — 
1.9658149(2) (2ll42^. For D < D t the behavior along the 
D — T critical line is governed by the Ising spin-1 jl fixed 
point at D — t — oo (at which vacancies are effectively 
forbidden) (HH]. 

Therefore, as regards the respective universality class 
of the (suppressed) phase transition, by fixing D < D t 
and taking ho ^ one is in fact crossing the critical line 
far away from its governing fixed point. Though this does 
not invalidate the conclusions of Ref. [l(| which concern 
the validity of the local-density approximation for cal- 
culating the internal energy and linear structure factor, 
one must expect strong crossover effects if attempting to 
evaluate critical indices there. 

In what follows, we kept D = 0, and p = 2 in Eq. (|17[). 

The "order parameter" associated to the trapping 
transition, i.e. the vacancy concentration, can be extrap- 
olated from TM data for 4 < L < 13 to be 1 - (S 2 ) = 
0.16120(1), in excellent agreement with the estimate from 
Table IV of Ref. [H namely p(Ising) = 0.1610(5). This 
is a reminder that one is working away from the critical 
point. Thus, for instance, in this case it is not possible 
to define the analogue of the magnetization exponent <p 
of SecIVl 

On the other hand, since for sufficiently small ho the 
associated length £ still diverges, at t = one can exam- 
ine the interplay between £ and, say, £l(0, £), and find the 
respective scaling behavior and (apparent) exponents. 
With this proviso, the arguments leading to Eq. (fl"5|) 
are still valid. Using 4 < L < 15, and 2 < £ < 9, we 
attempted data collapse of L/7t£l(0,£) against L/£ app , 
with the best superposition for 9 app = 0.700(5). The re- 
sults are shown in Fig. [7] Comparison of the respective 
numerical values shows that this scaling curve exhibits 
no superposition (except for the limiting constant value 
of 1/4 at x — >• 0) with that for the 5 = 1/2 Ising case, 
shown in Fig. [foJ 

In a similar vein, we investigated the site energies. 
Making now £/L = 1/V2, adjusting a sequence of data 
for 4 < L < 15 to e s = e° s + a'£ 9 '^, we get 0.10 < 
^ap P ^ 0.13. The discrepancy between 9 app and 8 app of 
the correlation-length scaling again reflects the fact that 
one is not close to an actual fixed point. 



VI. DISCUSSION AND CONCLUSIONS 

Thermodynamic systems with degrees of freedom sub- 
ject to a trapping field bring yet another length scale 
into play, namely, that of the trap size, £. The latter is 
a measure of how far from the center the trapping field 
reaches a given intensity, ho = £~ p , with p describing the 
steepness; p = 2 for the usual harmonic traps. There- 
fore, analyses of data extracted from systems of finite 
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Figure 7: (Color online) Scaled correlation-length data for BC 
model, p — 2, with varying L/l [see Eq. {T5}]. The horizontal 
axis is x = L/^" app , with 6> app = 0.700. The straight line 
through the origin is a guide to the eye. 

extent, L, (at least in one of its d dimensions), must 
incorporate the interplay between these two independent 
length scales. We have proposed an extension of the usual 
finite-size scaling theory to deal with these problems, and 
tested it on classical (Ising) systems, for which the scal- 
ing limit of large L can be easily achieved. The tests 
have been carried out on long strips of finite width L, 
subject to a trapping field translationally invariant along 
the "infinite" direction; the strip topology was that of a 
square lattice. 

For all field regimes considered, the characteristic ex- 
ponents of trap scaling vary continuously with p. This is 
analogous to what happens in the two-dimensional XY 
model, which has a line of fixed points, and exponents 
varying along the line [24j. There, the position on the 
line is related to a relevant variable which is marginal in 
the sense that it does not remove criticality, although it 
does affect exponents. Here we have the same scenario 
(but without, yet, a flow diagram to support it). 

In the low-field regime (i.e., for "shallow" traps, £ ^ L) 
one expects a linear-response scenario to be appropriate. 
For Ising systems we can therefore use this, together with 
exact conformal-invariance results, to produce analytical 
predictions to be compared with numerical results from 
TM methods, at the critical point of the otherwise infi- 
nite and untrapped system. We have established that, 
in accordance with our ansatz, the local magnetization 
displays the scaling behavior 

(a yi ) = h L^ +2 -rtf p ( yi /L), (18) 

to leading order in ho, where y\ is the distance from the 
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center along the finite strip direction; see Sec. Mil 

In the high- field regime, I < L, scaling bears the sig- 
nature of the p-dependent trap exponent 6, as previously 
discussed in Ref. |8| [see Eq. ©]. For p = 2 we have pro- 
vided independent checks of the exponents for the scal- 
ing at criticality of the central magnetization and of the 
energy density, in very good numerical agreement with 
Ref. |8|. Furthermore, with our TM calculations over a 
broad range of values of I and L (spin-1/2), we have es- 
tablished that the scaled correlation lengths behave as 
L/-7t£l(£) ~ g(L/£ ), where the function g(x) very ac- 
curately follows a form predicted by scaling arguments. 
Controlling for the quality of data collapse provides us 
with an estimate for the exponent 9 w 0.51 (for p = 2) 
which agrees very well with the respective scaling predic- 
tion, 6 — 16/31. We have also given an analytic treat- 
ment which describes how the finite-p scaling function 
evolves, in the p — > oo limit, towards a piecewise straight- 
line shape [see Eq. fll6|)]. 

As a further check of the theory, we have also consid- 
ered the high-field regime for a spin-1 Ising model in two 
situations. The first corresponds to an immediate ex- 
tension of the spin-1/2 case, in which the trapping field 
couples to the magnetization. Since the TMs for 5 = 1 
are larger than for S = 1/2, we were not able to consider 
linear lattice sizes as large as previously. Thus, the qual- 
ity of the data collapse and scaling fits was somewhat 
compromised; however, our results are broadly consis- 
tent with fits to the same exponents as before, namely, 
6 = 16/31 and <fi = Ofi/v = 2/31 (in standard notation 
of critical exponents) for p = 2. This picture therefore 
gives support to the universality of trapping exponents 



against spin magnitude. 

We have also examined the Blume-Capel (BC) model, 
with the trapping field coupling to £\ Sf (whose average 
gives the "particle density", in the lattice-gas language). 
This latter feature might be considered a more transpar- 
ent way to connect the spin language to that of particle 
trapping. However, we have seen that the fixed-point 
structure of the BC phase diagram puts one at a disad- 
vantage, concerning the estimation of scaling exponents 
for physically plausible values of the model's parameters. 
Notwithstanding this, we have been able to show that the 
competition between the associated scaling fields, L^ 1 
and gives rise to an effective scaling picture in qual- 
itative agreement with standard FSS ideas. 
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